The dihadron fragmentation function and its evolution 
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Dihadron fragmentation functions and their evolution are studied in the process of e'^e~ annihi- 
lation. Under the coUinear factorization approximation and facilitated by the cut-vertex technique, 
the two hadron inclusive cross section at leading order (LO) is shown to factorize into a short dis- 
tance parton cross section and a long distance dihadron fragmentation function. We provide the 
definition of such a dihadron fragmentation function in terms of parton matrix elements and derive 
its DGLAP evolution equation at leading log. The evolution equation for the non-singlet quark frag- 
mentation function is solved numerically with a simple ansatz for the initial condition and results 
are presented for cases of physical interest. 
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I. INTRODUCTION 

Lattice QCD calculations predict a phase transition 
from a hadronic gas to a quark gluon plasma (QGP) at 
very high energy densities in which quarks and gluons are 
no longer confined to the size of individual hadrons. To 
create such a dense and hot matter, heavy ions are ac- 
celerated to extremely high energies to collide with each 
other. If formed in such heavy-ion collisions, the QGP is 
rather short lived and hadronizes quickly into a plethora 
of mesons and baryons. Hence, the existence of such 
a state in the history of a given collision must be sur- 
mised through a variety of indirect probes. One of the 
most promising signatures has been that of jet quenching 
which leads to the suppression of high px particles 
emanating from such collisions. Such jet quenching phe- 
nomena have been among the most striking experimen- 
tal discoveries from the Relativistic Heavy Ion Collider 
(RHIC) at Brookhaven National Laboratory. Not only 
has the suppression of single inclusive high pT hadron 
spectra been observed Q, but also the disappearance of 
back-to-back correlations of high pT hadrons has been 
noted 0. Both phenomena and the observed azimuthal 
anisotropy of high px hadron spectra are qualitatively 
consistent with the picture of parton energy loss that 
leads to jet quenching. This is indicative of the forma- 
tion of a hot medium which is opaque to energetic partons 
and has a parton density about 30 times higher than in 
a cold heavy nucleus. 

In the investigation of jet suppression, correlations be- 
tween two high pt hadrons in azimuthal angle are used to 
study the change of jet structure While the back-to- 
back correlations are suppressed in central Au + Au col- 
lisions, indicating parton energy loss, the same-side cor- 
relations remain approximately the same as in p + p and 
d+Au collisions. Given the experimental kinematics, this 
is considered as an indication of parton hadronization 
outside the medium. However, since the same-side cor- 
relation corresponds to two-hadron distribution within a 
single jet, the observed phenomenon is highly nontriv- 



ial. To answer the question as to why a parton with a 
reduced energy would give the same two-hadron distri- 
bution, one has to take a closer look at the single and 
double hadron fragmentation functions and their modi- 
fication in medium. In this paper we take the first step 
by studying the dihadron fragmentation functions in the 
process of e^e" annihilations. 

Inclusive hadron production cross sections in e^e" col- 
lisions have turned out to be one of the many success- 
ful predictions of pcrturbative QCD 0, 0, 0- For re- 
actions at an energy scale much above h-qcD one can 
factorize the cross section into a short-distance par- 
ton cross section which is computable order by order 
as a series in as(Q^)] and a long-distance phenomeno- 
logical object (the single hadron inclusive fragmenta- 
tion function) which contains the non-perturbativc in- 
formation of parton hadronization 0|. These fragmen- 
tation functions can be defined in an operator formal- 
ism and hence are valid beyond the pcrturbative 
theory. They, however, cannot be calculated pcrtur- 
batively and have to be, instead, inferred from experi- 
ments. The definition of these functions affords them the 
mantle of being universal or process- independent. Once 
measured in one process, e.g. e+e^ annihilation, they 
can be applied to another, e.g. deep inelastic scatter- 
ing ov p + p collisions, and therein lies the predictive 
power. Another contribution of pQCD rests in the fact 
that once these functions are measured at a given en- 
ergy scale, they can be predicted for all other energy 
scales via the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) evolution equations [H 111 111 ■ 

The single inclusive fragmentation function D''^{z, Q^) 
can be interpreted as a multiplicity distribution for 
hadrons of type h with a longitudinal momentum fraction 
z that materialize from a fragmenting parton of flavour 
q. One can have, in principle, an n-hadron fragmen- 
tation function D^'^''*^' - ''" (zi, Z2, ...Zn,Q^) which counts 
the number of hadrons of type hi, h2, ...hn with mo- 
mentum fractions zi, Z2, ...Zn materializing from a frag- 
menting parton q. In this article, we will be con- 
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cerned with the double inclusive fragmentation function 
(zi, Z2, Q^) or the dihadron fragmentation func- 
tion. The operator definition of this function is not 
merely a trivial extension of the single hadron case; there 
are no straightforward sum rules connecting it to sin- 
gle inclusive fragmentation functions. Similar functions 
have been considered previously 0,0,0]. However, 
such formulations considered the double fragmentation 
function with a fixed angle 5 (or a fixed tranverse mo- 
mentum difference qt) between the two hadrons. For a 
large enough choice of 6 (or qx) the dominant contribu- 
tion to D^^''^'^{zi, Z2,5,Q^) was postulated to emanate 
from a process where the fragmenting parton had un- 
dergone a split into two partons which then fragmented 
independently. This formulation, however, requires S (or 
qt) to be large enough that the splitting process may be 
calculated in perturbation theory. In our formulation, 
the internal angle S (or relative momentum qx) will be 
explicitly integrated over. In this regard, our calcula- 
tions are similar in spirit to those of Ref. In that 
effort, general evolution equations for n-hadron fragmen- 
tation functions were motivated, and an algebraic solu- 
tion of the moments of the fragmentation functions was 
obtained. In none of the previous studies, however, has 
an attempt been made to generalize the operator defini- 
tion of the fragmentation function from one hadron to 
many hadrons. Evaluation of the n-hadron production 
cross section at leading order (LO) and leading twist al- 
lows one to analytically define such a function. Evalua- 
tion of corrections to the same cross section at leading log 
(LL) and leading twist allows one to derive the evolution 
equations for such functions. This is the objective of the 
current study. To the knowledge of the authors, no such 
calculation has hitherto been attempted. 

Our eventual goal is to derive the medium modification 
to the dihadron fragmentation function as the fragment- 
ing parton propagates through the medium. As shown 
in the case of single hadron fragmentation functions , 
the medium modification of the fragmentation function 
due to multiple scattering and induced gluon radiation 
resembles very much that of radiative corrections due to 
gluon bremsstrahlung in vacuum. Therefore, the study of 
DGLAP evolution of dihadron fragmentation functions in 
this paper can already provide us hints of what one may 
expect for medium modifications. 

The remaining sections are organized as follows: in 
Sec. II we present a general discussion of the double 
fragmentation function. We outline how such a function 
may be isolated in the expression for a double differential 
inclusive cross section and discuss the possible nature of 
its evolution equations. In Sec. Ill we begin with the 
S-matrix expression for the double differential cross sec- 
tion for the production of two hadrons in e'^e~ collisions 
at leading order (LO); we then factorize the expression 
at leading twist into the conventional hard part and the 
double fragmentation function. We also derive the rules 
for such an object in the cut- vertex technique of Mueller. 
In Sec. IV we write down the double differential cross 



section at next-to-leading order (NLO) and once again 
factorize it at leading twist into the conventional hard 
part and NLO correction to the fragmentation function. 
Using these we derive the DGLAP [ll IH |l2| evolution 
equations for the double fragmentation functions. In Sec. 
V we focus on the evolution equation for the non-singlet 
(NS) quark fragmentation function and solve its evolu- 
tion equation numerically. Finally in Sec. VI we discuss 
the results of our calculation and present our conclusions. 



II. THE PARTON MODEL 

In this section, we present a general discussion of the 
properties of a dihadron fragmentation function within 
a parton-model-like picture with coUinear factorization. 
However, all our assumptions will be demonstrated to 
hold explicitly in an operator formalism at leading log 
and leading twist in the subsequent sections. 

We consider the following semi-inclusive process 
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hi + h2 + X 



of e~^e~ annihilation. We consider two-jet events where 
both the identified hadrons hi and /12 emanate from the 
same jet. At leading order in the strong coupling this 
occurs from the conversion of the virtual photon into a 
back-to-back quark and antiquark pair which fragment 
into two jets of hadrons. 

In this scenario, at a large of the reaction, one 
may factorize the cross section of single inclusive hadron 
production as [Tsf 



da 

dz 



J2af[D'^{z) + D^^{z)] 



(1) 



with D'^{z) and D'i{z) as the single inclusive quark frag- 
mentation functions. The total cross section for the an- 
nihilation of an electron positron pair to a quark and an 
anti-quark, ctq* at leading order is 



is 



(2) 



Here, is the fractional charge of the quark in units of an 
electron charge, s is the square of center of mass energy 
of the e+e~ pair, and Nc = 3 is the number of colors in 
the fundamental representation of QCD. The fractional 
momentum z represents the the light-cone momentum 
fraction of the hadrons to the parent partons, i.e.. 



P/i • n 

z = , 

p • n 

where we use the notation of bold face letters rep- 
resenting four-vectors. The lightlike four-vector n = 
[n^, n~, n^] = [0, 1, 0] is taken conjugate to a given mo- 
mentum, as yet unspecified. 
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Similarly, one can expect to obtain the two-hadron in- 
clusive cross section from Eq. ^ by replacing the single 
inclusive fragmentation functions Z?^(zi) with the double 



inclusive functions ZJ^^"^ 



dzidz2 



We will not discuss cases in which each of the two hadrons 
emerges independently from each of the back-to-back 
quark and antiquark jets. 

In the parton model [T^l; at NLO with a single gluon 
radiative correction, the two hadron inclusive cross sec- 
tion in e'^e~ annihilation can be expressed as a convolu- 
tion of the fragmentation functions with the differential 
partonic cross sections (see Sec. 3.3 of Ref. 0|), 



ffa{zi,Z2, Q^) 



] dyD'^''''^{xi,X2)dxidx2 



da, 

dy 

d(Tq 

dy 
dag_ 
dy 

dOgq 

dy 

{hi,zi) {h2,Z2)]dxidx2 
^) dy [d^^^{x,)D^,Hx2)+ 



dyO'^^'^^ {xi, X2)dxidx2 
dyD^;^^^^{xi,X2)dxidx2 
dy [D'^^-{x,)D^^^{x2)+ 



{hi,zi) {h2,Z2)]dxidx2. (4) 

There are two distinct types of contributions in the 
above equation. The first one is determined by the 
two-hadron fragmentation functions of single partons 
that have the pQCD differential cross sections, daq/dy, 
daq/dy, dag/dy, and momentum fraction y. The second 
contribution corresponds to two, almost coUinear, par- 
tons (a gluon and a quark or antiquark with pQCD cross 
sections dagq/dy and dagq/dy) splitting from the same 
parent parton and then fragmenting independently into 
hadrons. In this case, the quark (antiquark) carries mo- 
mentum fraction y and the gluon carries 1 — y. One of 
the identified hadrons comes from each of these partons. 

The two hadrons hi and /12 have momentum fractions 
xi and X2 of their immediate parent parton which itself is 
endowed with a momemtum fraction y or 1 ~y. Relating 
the inside parton variables xi,X2,y to the outside hadron 
variables we obtain 

xi = zi/y,X2 = Z2/y 

for the first three terms and 

xi = zi/y,X2 = 22/(1 - y) 

for the last two terms. Using the pQCD partonic cross 
sections in the massive gluon scheme [ig] , one can obtain 



the double differential cross section for the production of 
two hadrons with momentum fractions Zi,Z2 as. 



d^a 
dzidz2 



q U + ^ I 



,{y)\og{Qyml)+aJ^^ (y) 



[D^q-^'^^{zi/y,Z2/y) + D\^'''^{zi/y,Z2/y)] 



a. 



^^P,^aM^og{Q'l7n'^) + aJI ^ {y) 

l-Zo 



X D'^'''^Hzi/y,Z2/y) + 



dy 



y(i - y) 



-P 



,(y) \og{Q^/m 



{D'^H^i/y)D'^H^2/ii-y)) 
Hzi/y)D'^H^2/ii-y))} 



+ 1^2. (5) 



Here the switch between the indices 1 ^ 2 is only meant 
for the last y integration. The P{y) functions are the 
regular splitting functions which contain both the real 
and virtual contributions and thus have no infrared di- 
vergences. The Piy) functions contain no contributions 
from virtual diagrams. The /(y)'s are the scheme de- 
pendent functions obtained in the massive gluon scheme, 
where nig is the fictitious gluon mass introduced to regu- 
late the collinear divergences. In the subsequent discus- 
sion we will focus on the leading log (LL) piece of the 
above expression. Thus one can drop the scheme depen- 
dent / functions. 

Note that up to this point we have simply retraced the 
sequence of steps in the evaluation of the radiative cor- 
rections to the single inclusive fragmentation functions 
in the parton model. What is new in the case of two 
hadron inclusive cross section is the contribution from 
the splitting into a quark and gluon followed by indepen- 
dent fragmentation. The log{Q^ /nig) in this contribution 
originates from an integration over the transverse mo- 
mentum q± of the quark and gluon emanating from the 
split. For very small values of q±, other higher order and 
nonperturbative processes become important that will in- 
validate the picture of independent fragmentation of the 
two partons. For q± » Aqcd, however, the higher 
order corrections will be suppressed and the quark and 
gluon will fragment incoherently in LL approximation 
(this was first pointed out in Ref. ^3)- ^ simple proof 
of this statement has been included in the Appendix. In 
this paper, we introduce a cut-off scale ^± that separates 
two regimes of two- parton fragmentation according to the 
value of q± : independent fragmentation for q± > fj,± and 
coherent fragmentation for q± < fi±. Unlike the factor- 
ization scales that we will discuss shortly, ^± is not intro- 
duced to rcnormalizc the fragmentation functions but to 
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define the pcrturbative (or non-perturbative) part of tfie 
difiadron fragmentation functions. It is quite analogous 
to the cone-size of jet definitions p^ . 

To simplify the discussion in this paper, we will con- 
centrate on the non-singlet fragmentation functions: 

D%f^{z,,Z2) = D'^^''^Hz,,Z2)-Dl^'''^{z,,Z2). (6) 
We also use the following convolution notations. 

f %A{zr/y,Z2/y)B{v) 

J a y 



functions. 




Note that the log in the second term can be separated 
into two pieces 



dy 



1/(1 - y) 



A{zi/y,Z2/{l-v))B{y). 



The bare fragmentation functions Dq^q^g{zi, Z2) in 
Eqs. Q and (jSj) are not as yet physical, measurable quan- 
tities and are scheme dependent, since the cross sections 
expressed in terms of them contain collinear divergences. 
One can however introduce renormalized fragmentation 
functions such that the double inclusive cross section can 
be factorizcd in the form of Eq. Q and is free of collinear 
divergences. 

Factoring out the e^e" annihilation cross section in 
Eq. jsj, we are left with the scale dependent physical 
fragmentation functions which should be free of collinear 
divergences. The "non-singlet" physical fragmentation 
functions are, 



\og{ml/^L\) = \og{my^^) + log(/iVMi), 



(9) 



with the first one containing the collinear divergence and 
the second piece defining the independent fragmentation 
of two collinear partons. In the case where /i < /ij^, we 
essentially have a function that depends on two scales: 
the collinear divergences that are extracted from the sec- 
ond term in Eq. ^ cannot be factorized out at a scale 
below However, we may chose to have > in 
which case the second piece in Eq. © is a finite constant 
that may be simply reabsorbed into the definition of the 
renormalized fragmentation function. With the factor- 
ization scale chosen above the physical scale ji^, we 
may now express Eq. ((SJ as 



D 



NS 



(10) 



D 



^1 ,/l2 

NS 



hi,h2 



(zi,Z2,Q ) = i^ATS 



(7) 



zi+22 



NS^g 



27r 
*2n" 



1-Z2 



log(QVMi) 
log(/ii/m2) 



-1 



Here the switch 1 ^ 2 is meant solely for the second and 
third term. Such exchange will be made implicit in the 
rest of this paper. We will also drop the limits of the 
convolutions in the notation for brevity. In the case of 
independent fragmentation, we have also split the expres- 
sions into a term that solely includes contributions with 
above the scale /.t^. The second piece contains con- 
tributions below /i^ and thus receives large corrections 
from higher order and non-perturbative processes. This 
piece will have to be absorbed into a redefinition of the 
bare dihadron fragmentation function. 

We now have to introduce the factorization scale fi 
and redefine the bare dihadron fragmentation function in 
terms of a renormalized one and the single fragmentation 



27r 



In what follows we will always compute in the region 
where fi > 

We may now substitute Eq. ifTTij) into Eq. and con- 
centrate on the leading order and leading log (LL) sector 
of the fragmentation functions {i.e. we only keep terms 
to order as{Q'^)log{Q^)) to get 



(11) 



+ D%^,D^^^-*i^ ^Pq^q, log(QVM') + - j • 

The above dihadron fragmentation function in the NLO 
still contains the bare single fragmentation functions in 
the contribution from two independent parton fragmen- 
tation. 
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Resum into ^ ,5 ) 



g' 



+ 



+ 



Or' 



+ 



Or' 



+ 




+ 




+ 




+ 



+ 




+ 



FIG. 1: Diagrams to be resummed iteratively to obtain the evolution of the fragmentation functions. 



If we consider higher order processes in which an ad- 
ditional gluon radiation takes place after the split but 
before independent fragmentation, as shown in Fig. ^ 
another coUinear divergence will arise. This is exactly 
the same as in the NLO correction to the single inclusive 
fragmentation functions. One has to introduce renormal- 
ized single liadron fragmentation functions at a factoriza- 
tion scale /.Ji 



D 



hi 
NS 



,(12) 



indicates the regular convolution notation, 
= / ^A{z/y)B{y). In addition, the renor- 



where, the (5 
i.e. A^B 

malized gluon fragmentation function is defined as. 



E 

q 



D 



27r 



Pa- 



Pn 



, \og{ml/fil 



,(13) 



where D'^jq represents the quark or antiquark fragmen- 
tation function and the sum includes all flavours. The 
factorization scale ^1 for the single fragmentation func- 
tions needs not be the same as the factorization scale for 
double fragmentation functions. 

With both the renormalized single and double hadron 
fragmentation functions, we obtain the leading log and 
NLO expressions of the double hadron fragmentation 
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functions 



(14) 



27r 



P,^,g log(QVM') + 



(Ml)* gp,^,,iog(gVM') + . 



In the above discussion we have set to be large such 
that there exists a hierarchy of scales ^^cd ^1 
Q^. The above factorization is valid in the regime in 
which the fragmentation functions are measured at a 
scale /i such that /i^ < /i << Q. In this limit we may 
also set ^1 = /i to define both single and double frag- 
mentation functions at the same scale. Note that the 
single fragmentation functions at the new scale ji differ 
from those at at higher order in as and thus the cor- 
rection due to these may be ignored in the leading log 
expressions for large enough . 

The remaining task in our parton model evaluation of 
the double hadron fragmentation function is to iterate 
the radiative process as shown in Fig. ^ Each of the 
circular blobs represents a fragmentation function at the 
scale /J,. Differentiating the series with respect to log((5^) 
followed by a reorganization of the various terms leads to 
the evolution equation: 



aiogQ2 



2^ 



Pn 



Pn 



(15) 



Within the framework of the parton model we can pic- 
ture the process as the free propagation of a parton fol- 
lowed by its fragmentation into hadrons of which two 
are identified. Fragmentation may be preceded by the 
radiation of multiple soft gluons (this is the top line in 
Fig-ffl) • Occasionally the parent parton undergoes a semi- 
hard split into two offspring partons which then propa- 
gate freely of each other and then fragment independently 
into hadrons and one hadron from each of these offspring 
is identified. Prior to their fragmentation, the offspring 
may radiate multiple soft gluons as well. 

The excercise in this section is based on the validity of 
our assumptions about the nature of the fragmentation 
process, especially on the validity of Eqs. H3I5|) . It was 
assumed that in progressing from single inclusive to dou- 
ble inclusive cross sections the parton model dynamics 
would remain the single leading behaviour and more im- 
portantly would lead to Eq. ((SJ). Such a proof exists for 
the single inclusive fragmentation functions that requires 
an exact operator definition of the single fragmentation 
functions. One can demonstrate both the factorized form 
of Eq. ^ and the evolution equations of the single inclu- 
sive fragmentation functions, exactly, as the leading log 



behaviour at large in an operator formalism (see Refs. 
H H 123 and Ref. Mounting such a proof for the 

dihadron fragmentation functions will require us to pro- 
vide a definition of the dihadron fragmentation function 
in the operator formalism. This will be the subject of the 
next section. We will extend the cut-vertex formalism of 
Mueller to dihadron fragmentation in this paper. 

The factorization of the NLO expressions at leading twist 
and leading log will be demonstrated in Sec. IV. 



III. THE SINGLE AND DOUBLE 
FRAGMENTION FUNCTIONS 

In the previous section, we made use of the parton 
model to motivate a double inclusive fragmentation 
function and assumed a factorized form as the leading 
behaviour of the two hadron inclusive cross sections. To 
prove the factorized behaviour we need to first obtain a 
consistent definition of the dihadron fragmentation func- 
tion. 

We begin with the matrix element for the electron 
positron annihilation into a given state of hadrons in the 
single photon approximation. 



-j:(ki+k2)-y 



(ki 



C(k2)7^^(ki). (16) 



In the above equation, J^iy) = J2q ^q''Pq{v)l'''^q{v) is the 
hadronic eletromagnetic current and ki , k2 are the mo- 
mentum four-vectors of the electron and positron. Here 
the sum over the number of colors in the fundamental 
representation of QCD is implied. Squaring the matrix 
element, summing over all final states of hadrons and av- 
eraging over all initial spins of hadrons, one obtains the 
total cross section for e"''e~ annihilation into hadrons, 



2s ^ J 2Ef{2TTY 

Shad 



(27r)4^(ki+k2-P5,„J 



Cf,Ao\J''mShad){Shad\r{o)\o) 



4(q2)2 M 
2sq4 4 



(17) 



where is the leptonic tensor and W'^'^ is the hadronic 
tensor. The four-momentum of the virtual photon is q 
ki + k2 = ((5,0,0,0) and the Mandelstam variable s = 
q2 = rpj^g ^^gj. g^^^^ includes both the sum of 

the complete set of states and the phase space integration 

n/e5.„.rf'P/W(2T)^- . 

One can evaluate the single inclusive cross section by 
summing over all possible hadronic final states that con- 
tain the identified hadron h. In the leading order and 
leading twist in a coUinear approximation, one can obtain 
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Eq. 1^. In a light-cone gauge (n ■ A = 0), the operator 
expression for the single inclusive frag mentation function 
at leading twist is obtained as [3, 12(1. , 



4 

X Tr 



(27r)4 
7 • n 



Ph ■ n 
p ■ n 



Ph ■ n 



Tq(q)ip,Ph) 



(18) 



where the the Dirac operators T^jj-^-) (p, p^) are given by 
(f,)„,^(p,Pft) = Jd^x J2 {0\MO)\Ph,Shad-l) 



X (p^,^-l|V;^(x)|0}e'P-^ (19) 



ifg)^AP^Ph) = Jd^x J2 {0\MO)\Ph,Shad-l) 

X (p,„^w-l|Va(x)|0)e'P-". (20) 

Here, the sums are taken over all physical final states 
of hadrons, which always contain, at least, the single 
hadronic state with momentum p^. In this case n is 
chosen such that its spatial components are antiparallel 
to the spatial components of the observed hadron. This 
implies n ph = = (p/J + |ph|)/2. In our choice of 
light-cone momenta p~ = p^ — pz- The gauge links re- 
quired to make this expression gauge invariant have been 
suppressed as they do not contribute to the leading twist 
fragmentation functions in light-cone gauge. 

The fragmentation functions can also be reexpressed 
in the cut-vertex technique of Mueller. These represent 
a powerful computational tool that may be used to calcu- 
late inclusive cross sections and the scale dependence of 
these functions in perturbation theory in a diagramatic 
languange. The Feynman diagrams illustrating the lead- 
ing order expressions for the single inclusive fragmen- 
tation functions are shown in Fig. |21 In this Feynman 
diagram the rule for the bare quark cut-vertex is 



7 • n / p,, • n 

6 \ Zh 

2p/i ■ n V P ■ n 



(21) 



The derivation of operator definitions for dihadron frag- 
metation functions and the extension of the cut-vertex 
technique to incorporate these functions is the focus of 
this section. 

We will concentrate on the two-jet events in 
e"'"e~ annihilation and are interested only in two hadron 
production off one single jet. The fate of the "back-side" 
jet will not be dwelled over here. We assume that the 
sum over all hadronic states in W^" can be simplified 
into two complete sets of states and that each overlaps 
independently with the quark and antiquark jet. This as- 
sumption neglects the interference between the two jets 
and is valid in leading log and leading twist. We also as- 
sume the duality between the complete hadronic states 



Vh 




FIG. 2: cut- vertex for quark fragmention function at LO. 



and partonic states. Thus the sum over hadronic states 
in the "back-side" jet will be replaced by partonic states. 
Under this assumption, wc evaluate the hadronic current 
operator, order by order in a^, by expanding the QCD 
interaction Hamiltonian in the interaction picture. In the 
leading order (LO), one obtains. 



E E 



d^pid^P2 



AEiE2{2TTf 
(27r)4,54(q-pi-p2-k-^ p/) 



d^k 

^ ' feShad-2 
X (fc,pi,p2, Suad - 2|i^",(0)7''V,(0)|0). 



(22) 



In the above equation. Shad is a complete set of hadronic 
states. In the remaining discussion we will drop the sub- 
script {had)- We have extracted two particular hadronic 
state sums, labeled as Pi,P2, from the full sum over states 
S. The hadronic tensor may be represented by the Feyn- 
man diagram in Fig. |31 

On Fourier decomposition of one of the quark or an- 
tiquark operators, followed by a sum over all spins of 
the outgoing antiquark (quark) state k, we obtain the 
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where the parton-double-hadron overlap matrices are 
given similarly as for single fragmentation functions as 



af) 



S-2 



{Q\mx)\pw2S - 2){p,p2S - 2|V^^"(0)|0) (25) 



and 



Tq{p]Pl,P2) 



S-2 



{0\i;^ix)\piP2S - 2){p,p2S - 2|V;,^(0)|0), (26) 

where the two quark field operators have the same color 
index and average over colors is explicitly implied. 

Up to this point the derivation of the factorized dou- 
ble inclusive fragmentation function has followed a path 
not entirely dissimilar to that of a single fragmentation 
function. At this point one may introduce the light-cone 
variables and their ratios, 



FIG. 3: The leading order Feynman diagram contributing to 
the double inclusive fragmentation function. 



Zl = Pt/P^ , Z2 = pt/p'^ 



hadronic tensor as. 



4-/^7 AE,E2{2tt 



9, s-2 



(23) 



(27r)' 



{(0|^,(0)|pi,p2, S - 2)7^' ^7"(P1,P2, S - 2|^g(0)|0)} 

+[rmm\pi.P2. s - 2){p^,p2,s ~ 21^5(0)10)7'' ^} 



One may rewrite the 5— function as a four-space inte- 
gration of an exponent exp(— ip • x) that in turn can 
be used to transform the quark wavefunction operator 
■0(0) — > ip{x). A shift in the d'^k integration can be per- 
formed i.e. k = q — p. Summing over the spins (and 
colors) of all the final (parton) states we obtain the gen- 
eral form of the hadronic tensor as 



(Ppi(Pp2 f d^p 



2^5+((q-p)2) 



Tr 



AEiE2i2T:)^ J (27r)4 
{fq{p;p,,p2h''U- ^h"} 



(24) 



+ Tr {fgip; p,, p2)r U~ ^h''} 



and z = zi + Z2 = {pi + pt)/p^ = Ph/P~^- Where ph 
is the total momentum of the pair of hadrons that are 
identified. We now take the coUinear approximation that 
the hadron momenta pi,p2,Ph are almost collinear with 
respect to the quark (antiquark) momentum at high en- 
ergies and are thus dominated by their + components 
for light-cone vector n chosen in the direction of the out- 
going quark or antiquark. Essentially, p^ » Pi,pi± 
and the same is true for p2 and ph. The overlap ma- 
trices Tq, Tq have a Dirac matrix structure and hence 
can be decomposed in a basis of products of 7 matri- 
ces (1, 7'', cr^'', 7^7^, 7^). The only term of this basis to 
survive is 7^: even combinations are set to zero in the 
trace and combinations containing 7^ vanish under a spin 
sum. As Tq, Tq are scalers and only depend on the three 
almost collinear momenta, we obtain at leading twist the 
following decomposition of the overlap matrices: 



Tq{p;pi,P2) = ^rg(p;pi,p2) or 



Tq{p]Pl,P2) 



Tr[rlfg{p;pi,P2)] 
2n ■ pk 



(27) 



where Tq{p;pi,p2) is a scalar function. The entire Dirac 
structure has been extracted into the 7 matrix. 

With the aid of (5— functions, we can introduce the 
definition of the fractional momenta zi and Z2 into the 
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hadronic tensor: 

= [ dzidz2d{l ~ zi - Z2) 
Jo 



.(28) 



Tr 



AEiE2{2T:f 

( Ph al , ,y\ 



Pt 



Pt 



j^Tgip; ,pi,P2)5{zi - -^)S{z2 ~ ^) 



d^p 



Tqip; ,Pl,P2) 



p""" p^ 



In collinear approximation, we expand the hard part 

S+iiq-pf). (29) 



H 



lliu 



Tr 



in the transverse momentum of the hadrons phi_ and take 
only the leading term T-L'^^{p'^) ^'H'^^{p'^ = p\l z)^ .... 
This approximation allows us to factor out the hard 
part from the d^p integral. This is the first step in the 
factorization of the double hadron inclusive cross sec- 
tion. Based on the collinear approximation we have also 
dropped the term or v\l z^ from the argument of the 
(5+ function (i.e. — 2q • p >> p^). Given the light- 
cone structure of the four- vector p a further simplication 
of the argument of the 5— fmiction may be obtained: 



reorganized as follows, 



LZ2 J (Ztt)-* Zi Z2 

1 r/P^ Pt\ f d'^P +2„, + pt . 



Z1Z2 ' Zi Z2 ' J (27r)4 



1 x/Pi Ph ^Pi \ f P 



Z1Z2 Zl 



Z2 



'Sipt - ^ 



Zl + Z2 



Zl + Z2 J (27r) 



(2^)4 



,P^ Sip+ - 



Zl 

P± 

Zl 



z z J (,^7rj z 

..pl5{z2pt^zipt)j ^Jiz^^l)...^ 



(32) 



where z = zi-\- Z2. Substitution of these expressions into 
the hadronic tensor W^^'^ followed by a substitution of 
W^^'^ back into Eq. H17() leads to the following expression 
for the total cross section for e+e~ annihilation. 



(T = e^'A^c j dzidz20{\ — Zl — Z2) 



X / dpidp22TT—5+{ph - zQ/2)p+d{z2pt - zip+) 



^Siz-Ei. 



Tqip;Pl,P2) +Tq{p]pi,P2) 



(27r)- ■ p 

where the external (5— function has been used to set 



(33) 



q' 2q • p 


= q^q 
= Q'- 


(Q - 


Q' - QP+ 


= Q^- 


9^, 



Qp' 



- q± ■p±) 

as q± = 



2q 



Ph_ 

z 



(30) 



With the above simplifications and reorganization of 
arguments in the two internal S functions, the hadronic 
tensor can be written as 



Nc I dzidz2c 
'0 



Zl 



Z2) j dpidp2 



(31) 



A 



2Tr6+{Q^ ~2Qi-^)n^"' 

z 

2 / d'^P 

9 



{27rr 



Tqip;,Pi,P2) +Tg(j);,pi,P2] 



p 



Z1Z2 



S{p+ P^)S{p+ - ^), 

Zl Z2 



where we have introduced the shorthand notation dp = 
d^p/ {2TiY2Ep. The variable p^ is ovcrdetermined and 
thus one of the 5 functions acts really on the integrations 
external to d'^p. These may be extracted and further 



2Q 



5+{Ph-zQ/2)W-{Q). 



The dependence of the hard part on hadronic variables 
is replaced with the appropriate partonic variables. Dif- 
ferentiating the above equation with respect to zi , Z2 
leads to the double differential cross section outlined in 
Eq. (jSJ). Before the extraction of the double inclusive 
fragmentation function, some simplification of the above 
equation is in order. In contrast to the definition of the 
single fragmentation function, there are double hadronic 
integrals d^pi,d^p2 and two sets of iS— functions as op- 
posed to one. The cautious reader will note that the over- 
lap matrices Tq{p\pi,p2),Tq{p\ ,pi,p2) [see Eqs. H25l26|l ] 
arc dimensionally different from the overlap matrices 
in the definition of the single inclusive functions (see 
Ref. [111). 

To simplify, we begin with a variable trans- 
formation. One essentially changes from 
the set [pi^,piy,pi^^p2^,P2y,P2z] to the set 
[pi,P2,qi.,0cm,4>cni,4>\] as illustrated in Fig. |3| This 
choice is not entirely arbitrary. The discussion of the 
NLO in the previous section required us to incorporate 
transverse momenta q± up to a semihard scale into 
the bare fragmentation function. This particular variable 
transformation allows us to isolate the q±^ integration. 

The new vector pf^ ~ P1+P2 has the three components 
of mostly massless four- vectors. The requirement that 
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FIG. 4: 



The variable transform 

i/.J'izi-P2j;,P2y,P2,] to [pi,P2,q±,0cm,(t>c,n,(t)±] 



from 



= Pi +P2 is trivially fulfilled. The new "body- fixed" 
variable q± quantifies the component of pi — p2 that lies 
on the plane perpendicular to ph and (j>± is the azimuthal 
angle of q± on this plane. The angles 9cm and (ficm quan- 
tify the direction of ph with respect to the e"'"e~beam 
direction as the z axis. The Jacobian for this tranforma- 
tion is simply 



at leading twist. With these new variables one may eas- 
ily relate the partonic variable in the hard part with 
the hadronic variable ph as p/z, i.e. the sum of the 3- 
momenta of the detected hadrons is collinear with the 
3-niomenta of the fragmenting quark or antiquark. It 
may be demonstrated that the corrections to this state- 
ment contribute at higher twist. It should be pointed 
out, in passing, that this is a more accurate statement 
than the assumption of coUinearity between the leading 
hadron and the fragmenting quark (antiquark) in the case 
of single fragmentation, as in most cases a dominant part 
of the jet's momenta is contained in the momenta of the 
leading hadrons. 

With these new variables we can evaluate the inner 
product of the leptonic tensor and the hard part of the 
hadronic tensor, 



and obtain the double differential cross section as 
da ^.^e^elNc f dpidp2 



dzidz2 



e{i~z,-z2)Y, 



d cos 6, 

X (1 + COS^ 0cm) 



2s 



i2n)HpiP2 



dq±d4)±^{pi +P2f 







fp_i 




(^'^ 4) 






Z1Z2 




Z2 ) 



d'^p 

(2^ 



5 z 



Pt 



Tq{p;pi,P2) + Tqip;pi,p2) 



(35) 



Assuming that the overlap matrices, Tg(j);pi,p2) and 
Tq{p;pi,p2), are independent of the angles, 9cm, 4>cm and 
(j)±, one can carry out the integrations over these vari- 
ables. It may once again be stipulated that the following 
factorization is being performed and expected to succeed 
only at high energies and momenta (large limit). In 
this limit we note 



P 



• P 



i-o(4 



Hence, in the collinear limit p,p^ >> q±, one may sub- 
stitute p p^ in the entire integrand. Regardless of the 
presence of S functions it may be demonstrated that the 
correction to this approximation is suppressed by at least 
a power of Q^. The two remaining i5— functions external 
to the d'^p integration, i.e. 



may be used to evaluate the pi and p2 integrals. Factor- 
ing out the LO total e^e^ annihilation cross section CTq^, 
we obtain the factorized double differential cross section 
in a form similar to the structure of Eq. Q as. 



da 



dzidz2 



^ \^ qq f q±dq± f 
4-''° i 4(27r)2 4ziZ2 7 



d^p 
(2^ 



(36) 



Tq{P;Pl,P2) + Tq{p]Pl,P2) 



Slz-P^ 

p+ 



We thus arrive at the definition of the leading order 
double inclusive fragmentation function as 



dqj 



8(27r)2 AziZ2 J {2tt 



d^p 



(37) 



X Tq{p; Zip, Z2p)5 [ z 



Pt 



P^ 

In cut-vertex notation, the dihadron fragmentation func- 
tion may also be expressed by the following equation 



A-^''" -4g'*(l+C0s2 9cm) 



(34) 



D'i'''Kzi,Z2) 



4:ZlZ2 



Tqizi,Z2), 
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where T{zi,Z2) is given by the diagram in Fig. |S| Note 
that the bare cut-vertex has undergone no change as com- 
pared to the single hadron fragmentation function, ex- 
cept it takes as input the sum of the fractional momenta 
z = zi + Z2. The soft hadronic sector is slightly modified 
by the exclusion of two hadronic momenta (instead of 
one) and the integration of the transverse momenta qj_ . 




FIG. 5: The cut-vertex representation of the dihadron frag- 
mentation function. 

Note that the definition of the dihadron fragmentation 
functions in their factorizcd form as in Eq. ^ seems 
to depend on our choice of variable transformation (sec. 
Fig. El- It nia-y indeed be possible to use a different vari- 
able transformation and obtain a similar factorization. 
The sole constraint on the choice of transformation to be 
used is based on the factorization of the NLO expressions 
into a form similar to that of Eq. ||SJ), with the same def- 
inition of the double fragmentation functions. We will 
demonstrate this in the next section. A second motiva- 
tion for this choice of variables is the ability to isolate 
the transverse momentum q±, which is integrated over. 
In the operator definition of the dihadron fragmentation 
functions. q± is generated non-pcrturbativcly. So wc will 
call it intrinsic. The upper limit of integration of this 
intrinsic q± has not been specified in Eq. 1)3 7|l . From our 
discussion in the previous section of the parton model 
and of the NLO processes in the next section, the up- 
per limit may be set as We will assume that hadron 
pairs with the relative transverse momentum q± > fij_ 
are generated only perturbatively. 



IV. CROSS SECTION AT NLO AND DGLAP 
EVOLUTION 

With the definition of the dihadron fragmentation 
functions in the operator formalism, which is shown to 
factorize from the hard parton cross section in the LO, we 
are ready to study the DGLAP evolution of the dihadron 
fragmentation functions by computing the double inclu- 
sive cross section at next to leading order (NLO). The 
calculation will also justify the factorizcd form of Eq. 
(0) and our LO definition of the dihadron fragmentation 
functions. 

The NLO matrix element of e"^e~annihilation process 
can be obtained from the perturbative expansion of the 
S'-matrix with an interaction Hamiltonian which includes 
an interaction potential corresponding to a quark of color 
j interacting with an antiquark of color i (or visa versa) 
and a gluon of color a: 

Us = igtlj^^Y^,A';>- (38) 

The tfj's represent the GcU-Mann matrices. 

The DGLAP evolution arises from coUinear gluon 
bremstrahlung in the final state of two jet events. As 
in the case at LO in the previous section, we can again 
neglect the interference between the two jets in either 
the partonic or the hadronic level in the leading log and 
leading twist approximation. Therefore, the sum over 
all hadronic states is decomposed into the sum over two 
complete sets of hadronic states, each overlapping with 
one of the opposite-moving jets. Invoking the parton- 
hadron duality, the sum over the hadronic states in the 
"back-side" jet can be replaced by a sum over partonic 
states. Therefore, partonic processes within the "back- 
side" jet will not contribute to the evolution of the frag- 
mentation function in the opposite side which has been 
defined within a rather strict coUincar approximation. As 
in the case of single fragmentation function, wc will also 
assume, in addition, that there is no interference between 
the fragmentation of the leading parton and the radiated 
gluon that has a minimum transverse momentum set by 
the factorization scale /i. 

Hence, in NLO, the sum over all states may be ex- 
pressed as 

1^) = |fc,(,-)> X |5-2,pip2). (39) 

Depending on how different operators are contracted 
with the outgoing hadronic state \S — 2), we may identify 
three different cases: 

|^-2,piP2) = |^-2,piP2) X \{p~l)g) 

+ |^-2,PlP2) X \{p-l)qiq)) 

-f 1^-1, Pi) X |S'-1,P2). (40) 

These cases differ in the partonic operator that contracts 
with the hadronic state. In the above, the parton with 
momentum k which proceeds in a direction opposite to 
that of the identified hadrons alternates between a quark 
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and an antiquark. In the rest of this paper, we will al- 
ways take it as an antiquark in order to focus on the quark 
fragmentation function. The case for the antiquark frag- 
mentation function will be formally identical to that of 
the quark. 

In the following subsections, we will evaluate contri- 
butions from these different cases in detail. Roughly 
speaking, the first hue of Eq. (|^ represents contribu- 
tions where the fragmenting quark undergoes a split into 
a quark and a gluon and the two identified hadrons em- 
anate from the quark offspring. The second contribution 
represents the case where both identified hadrons em- 
anate from the gluon. The last contribution represents 
the case where one hadron emanates from each of the 
quark and gluon offspring. 



A. NLO contribution from quark fragmentation 

Proceeding with the evaluation of the double differ- 
ential cross section at next-to-leading order, the focus 
in this subsection will be on isolating contributions to 
the cross section that contain an explicit expression for 
a dihadron fragmentation function of a quark. The out- 
state in this subsection will solely be restricted to the 
first line of Eq. H40() . The instate is simply that of an 
incoming e^e~pair. Insertion of the interaction operator 
density T\HeJre-'y'Hqq-y'Hqqg\, foUowcd by a contraction 
of the outgoing antiquark operator with the states \kq) 
and gluon operator with \ {p — l)g), leads to the following 
matrix element: 



X (piP25-2|v;,(0)|0). 



-I- le 
(q - 1)2 + le 



(q - k)2 + ie 
X (27r)454(q-k-p-l) 



(41) 



In the above equation there are two terms with differ- 
ernt momentum dependences within the curly brackets. 
The reader will readily note that the second term is the 
Feynman rule for the process indicated in the upper panel 
of Fig. ini while the other term consists of the Feynman 
diagram where the gluon is radiated from the antiquark 
line, as shown in the lower panel of Fig. 

In computing the NLO cross section, one may once 
again factorize the cross section into a leptonic and a 
hadronic piece [see Eq. (|17|l ]. Summing over all final 
states of the outgoing antiquark, gluon and hadrons from 
the fragmenting quark (besides hi and /12) followed by 
an incorporation of minor simplifications, the hadronic 



k 




FIG. 6: The leading log contribution to the NLO modification 
of the quark fragmentation function. 



tensor may be expressed as 



W 



d'^p 



[2tt)HEiE2 J (27r)4 J (27r)4 ^"^^^ 
X 2^<5+((p - l)2)27r(5+((q - Y>f)g'' N.Cpdp^ip - I) 




(l-q)2 + k 



In the above equation the color factor NcCp comes from 
the factor Tr[^"^''](5°^ For brevity, we have omitted in the 
above the sum over quark flavors weighted with fractional 
charge, e^. The sum over polarizations of the gluon 
leads to the factor dp^. In light-cone gauge dp^ is given 



dpaip - = 9 per - 



(p - l)pna + (p - l)anp 
(p - 1) . n 



(43) 



The overlap matrix clement Tq(l;pi,p2) is the same as 
defined in the previous section, with the final fragmenting 
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quark momentum reduced to 1. The remaining portion of 
Eq. H42|) is also easy to trace. The two S functions essen- 
tially stipulate that the gluon and antiquark be released 
onshell. The first set of terms in the curly brackets repre- 
sents the process indicated by upper panel of Fig. The 
second set of terms indicates the case where the gluon is 
emitted from the outgoing antiquark. with all other fea- 
tures remaining unchanged, e.g., the fragmenting quark 
has momentum 1, the gluon still has a momentum of p — 1. 

Within the coUinear approximation, we again assume 
that ,pj[ » l~ ,pj^ ,l±,q±. As a result, the overlap 
matrix element may be factorized via the following ap- 
proximation: 



•' S-2 



X (0|^,(a:)|pi,p2,5-2)(pi,p2,5-2|V',(0)|0) 



^r,(l,pi,p2) 



(44) 



Within the coUinear approximation, we assume that the 
momentum of the final fragmenting quark is coUinear 
with that of the hadrons emerging from the parton frag- 
mentation. We thus define a new fractional momentum, 
z' = pj^/l^. This allows the partonic four- momentum 
vector 1 to be replaced with the hadronic four- vector: 
i.e. 1 = p/i/z'. 

As a result the factor (1 — q)^ is approximated as: 



(1 - q)2 ~ g2 _ 2QI+ ~ Q2 _ 2Qp+Jz'. 



(45) 



Unlike the case of the LO process, the transverse mo- 
mentum of the quark which emanates from the electro- 
magnetic vertex is non- vanishing, p_L ^ (this is also the 
transverse momentum carried by the gluon). As a result, 
the negative longitudinal momentum is constrained 
by one of the (5+ functions as. 



expression for the hadronic tensor, we obtain: 

X g''N,CF5{z^-p+lp+)5{z2 - pt /p+)dp,{p - I) 
X 2TTS+{2p-{p+ - 1+) -pI)2tt5+{Q^ - 2Qp+) 

7" , Y{y- ^)7"1 



X TV 



pll+/{p+-l+) Q^-2QV 



^^2 \pll+ 



7" ny- 

/ip+ - 1+) 02 _ 2QI+ 



- i=Pfi/^' 



(gj^(-'-f^)r,(l;pr,P2). 



A careful study of Eq. (|47|l reveals that the leading log 
contributions are dominated by the region where p± — > 0. 
As a result, the part of VF"" which represents the square 
of the process in the lower panel of Fig. El has no lead- 
ing log contribution. The leading log contribution comes 
from the square of the first term corresponding to the 
square of the process in the upper panel of Fig. El In the 
light-cone gauge, the interference terms have no contri- 
bution at leading log as in the case of single fragmenta- 
tion functions (see chapter 3 of Ref. 0|). This can be 
demonstrated to hold at leading twist merely by com- 
pleting the trace as indicated in Eq. (|47|l and extracting 
the p\ dependence from the numerators. 

The remaining factorization into the hard and soft 
piece proceeds as in the LO case, leading to the hadronic 
tensor at NLO at leading log and leading twist. 



dzidz2dz' 



d^Pid^P2 



d'^p 



{2TTfAEiE2 J (27r)4 
X g^N,CF5{zi-p+/p+)5{z2~pt/p'^) 
2t:5+ (p- 



X Tr 



2{p+ ~ 1+) 



2ttS+ [Q^ - 

7" y{z'/2)Y i>Y 

[pi? + /(p+ - Z + )]2 



1=1 



dH 



(,^),^(.'-f^)r,(l;p„P2). 



(p-p„/z')' - 2p+p--pi-2p-p+/z' = 

„2 



^ P = 



Pa 



2{p+-pI/z') 



(46) 



The trace of the Dirac matrices can be carried out by 
brute force. The leading-twist part can be obtained in 
a more straightforward way by rewriting the matrices 
within the trace symbolically as 



The hadron fractional forward light-cone momentum is 
stUl defined as z = p'^/p'^ as in the LO case. At leading 
twist one may replace all occurences of ph with z'l. The 
corrections to this approximation are down by powers 
of Q^. Incorporating the above approximations into the 



The leading twist portion of the matrix C, may be writ- 
ten as 

C=^C, 
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with the factor C expressed as a trace: 



C = Tr 



(49) 



i=Ph/2' 



where 



p' = p\}^ / {p^ ^ ^^)- This procedure is very sim- 



ilar to the one used to construct the colhnear approx- 
imation to the overlap matrices T'q(p;pi,P2)- One can 
complete the above trace and obtain the regular splitting 
function (1 -f y^)/(l — y), where y — z' j z. It represents 
the probability for the radiation of a giuon from a quark 
prior to its fragmentation. 

Using the above approximation one may extract the 
same hard part as in Eq. (jSHJ, 



Siiq-pH/zf) 



H"''<5((q-pjz)2). 



(50) 



Extracting this hard part and comparing with 
Eqs. (|28I31|I and the resulting cut-vertex diagram 
in Fig. Owe note that the soft part of Eq. H48|) begins to 
display a structure as illustrated in Fig. [7| 

To complete the calculation and obtain a factorized 
form of the NLO contribution to the double inclusive 
cross section, the integration over the tranverse momen- 
tum of the identified hadrons will have to be factored 
into the fragmentation function. Similarly as in the case 
of LO calculation, one has to transform the basis of the 
momentum integrations of the two identified hadrons to 
the basis indicated in Fig. ^ 

In the steps leading to Eq. H48|l . the approximation 
that the momenta of the hadrons and the fragmenting 
parton are colhnear has been made. Particular among 
these are the approximations that Z+ >> These 
essentially indicate that the invariant mass of the final 
fragmenting quark is negligible compared to its forward 
light-cone momentum. This is identical to the approxi- 
mation made on the momentum p in the LO calculation 
(see Fig. Since the NLO process has a colhnear di- 
vergence when p"^ 0, we will only consider the leading 
log contribution. 

We thus obtain the following factorized form for the 
hadronic tensor: 



dzidz2 



dpidp2 



d cos 0, 



X g^N.C, f% I ^E^|j^2.H"''5((p-l^) 



y 



Cp+H{pt - ZiP+)S{p+ - Z2P+) 



2y 



dqj 



(^^(^'-^)T,(l;P.P.). (51) 



Further factorization of Ti."'^ from the d*p integration 
leads to the factorized form for the NLO correction to the 
dihadron fragmentation function as illustrated in Fig.Qin 



pj+pt 



Pi 



P2 



4(2^) 



2p- 





FIG. 7: Next-to-Leading order cut-vertex for quark fragmen- 
tation and its factorization. 



terms of cut vertices. Inserting the hadronic tensor into 
the expression for the double differential cross section 
with minor simplifications, we obtain the factorized NLO 
contribution to the double inclusive cross section, 



dzidz2 



47r a'^e^Nc 

^ 1 02 



(52) 



2tt 



dpi 

2 



dql 



dy^ 1- 



A{z^/y){z2ly) J J (2^) 



y 
dH 



Ph. 

1+ 



Tq{l;pi,p2). 



The last line of the above equation may be easily iden- 
tified as D{zi/y, Z2/y, ^jl'^); the dihadron fragmentation 
functions scaled up by the momentum fraction y, carried 
by the quark emanating from the split. Physically this 
represents the contribution to the fragmentation func- 
tions at a higher order brought about by gluon radiation 
carrying away with it a momentum fraction I — y. By 
itself this process displays both an infrared divergence as 
y ~* 1 and a colhnear divergence as p± — *■ 0. The in- 
frared divergence will be canceled by the virtual diagram 
contribution as shown in the next section. The colhnear 
divergent part will be combined with the colhnear diver- 
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gent part of gluon fragmentation and absorbed into the 
renormalized fragmentation function. 



B. NLO contribution from gluon fragmentation 



We now proceed with the contribution from gluon frag- 
mentation in the NLO processes. Essentially the out- 
going hadronic state is replaced with the second line of 
Eq. H40() . Insertion of the interaction operator density 
followed by the contraction of the quark and antiquark 
operator with the out-going states \kq) and \{p — l)q) , 
leads to the following matrix element: 



(53) 



X (pip25-2|A^(0)|0)r(p-/)i;4^ 



(k + l)2+ze 



!.t;"(fc)(27r)'^(5(q-k-p-l). 



k 




FIG. 8: The leading log gluon fragmentation contribution to 
the NLO modification of the quark fragmentation function. 



As in the previous subsection, the approximation of 
very high energies is made, allowing the isolation of the 
leading twist and leading log contribution of the corre- 
sponding hadronic tensor. This is obtained as 



dpidp2dcos 9cmd4>c 



d'^p dH 

(27r)4 (27r)4 



27r5((q - p)')2^5((p - \f)g^ N,CfTS:Pi.P2) 



X dp^Tr 



12-26 



(54) 



where the gluon overlap matrix element Tg{l;pi,p2) is 
defined as 



dql 

X (pi,p2,5-2|4(0)|0) 



^^x (0|A;i(.T)bi,p2,5-2) 
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(55) 



ential cross section: 



dzidz2 



y-^ 47r a^e^qNc 



(56) 



27r 



dp2 



!_ I dy 1 + (1 - y)^ 

PI Jz y V 



{z/yf 



dq 



dH 



2{z,/y){z2/y) J J i^^) 

x<5(^^-^)t,(/;pi,P2). 



It may come as no surprise that the above equation 
may also be derived from a set of Feynman rules involving 
cut- vertices. The cut- vertex diagrams are illustrated in 
Fig. 1^1 The rules are indicated in the figure. As a result 
of this computation we may now present the cut-vertex 
expression for the gluon dihadron fragmentation function 
(indicated in the lower right hand corner of Fig. inj: 



A diagrammatical representation of this gluon frag- 
mentation process can be illustrated as Fig. |H1 The pro- 
cedure leading to the extraction of the leading log and 
leading twist is rather similar to the case for the single 
fragmentation function and to the case of NLO process 
of quark fragmentation in the last subsection. Contract- 
ing the hadronic tensor with its leptonic counterpart we 
obtain the gluon contribution to the NLO double diffcr- 



Dg{Zl,Z2) 



X 6 I z 



dqj 



dH 



2ziZ2 J 8(27r)2 J (27r) 



Ph 



^+ / Tg{l;pi,p2), (57) 



where the factor Tg{l;pi,p2) is given in Eq. (|55|) . 
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FIG. 9: Next-to-Leading order cut- vertex representation of 
the gluon fragmentation contribution. 



C. NLO contribution from quark and gluon single 
fragmentation 

In the previous subsections, we have evaluated two 
separate contributions to the NLO double fragmentation 
functions. In both of these we have assumed that the rel- 
ative transverse momentum of the two detected hadrons 
is intrinsic and is limited by a scale In the defini- 
tion of the dihadron fragmentation functions, the hadrons 
are detected with given fractions of forward longitudinal 
momentum but the transverse momenta are integrated 
over. Thus, all allowed transverse momenta between the 
detected hadrons must be included. In this paper, we 
will assume that all hadron pairs with relative transverse 
momentum larger than ji > fij_ are generated perturba- 
tively. In the next-to-leading order, such hadron pairs 
can be produced from the independent fragmentation of 
the quark and gluon after their split, as illustrated in 
Fig. ^1 Such a scenario has been considered in Refs. 
fl3, 14. 15| where the double inclusive cross section with 
two detected hadrons in e+e~collisions with a fixed trans- 
verse momemtum between them was computed. The au- 
thors argued that in the case that the transverse momen- 
tum lies in a semihard region Aqcd << Qi^ << Q, the 
dominant contribution to the cross section comes from 



the process where the fragmenting parton undergoes a 
semihard split into two independent partons which then 
fragment independently. 

Under the condition that << Q these hadrons cair 
still be considered to belong to the same jet. Moreover, 
when Aqcd << q±, the higher order contributions from 
multiple gluon vertex corrections to the semihard vertex 
are non-leading and thus the fragmentation of the two 
partons emanating from the split may be considered as 
independent. In Ref. the authors demonstrated that 
the Sudakov double logarithms from the higher order ver- 
tex corrections are absent in the region Aqcd << q±- In 
the interest of completeness we will repeat this derivation 
in a slightly different language in the Appendix. 

We will evaluate the leading log and leading twist con- 
tribution in which the hadron pair comes from the in- 
dependent fragmentation of the quark and gluon after a 
semi- hard split. We start again with the matrix element 
for this process. 



Qvfl - a 



X (pi, ^1 - 1|V;,(0)|0)(P2, ^2 - l|^;(0)|0) 



k)2 



(PSi 



It 



v(k) 



X (2^)^<5*(q-p5, -PS. -k) 



(58) 



The second term inside the square bracket corresponds 
to gluon emission from the antiquark in the quark direc- 
tion and will not contribute in the leading log approxi- 
mation. The out-state in this case is chosen to be the 
last line of Eq. H40|l . In this case, the sum over all final 
states of hadrons from the initial quark has been bro- 
ken into two identical complete sets. In each of the sets, 
S\ and 5*2, a single hadron will be identified. Unlike 
in the previous two subsections, where one of the quark 
or gluon operators is contracted with a partonic state; 
both the quark and the gluon operators will be contracted 
with hadronic states in this case. The cross section con- 
structed from this matrix element contains two separate 
sums over hadronic states; one of which has an overlap 
with a quark state, the other with a gluon state. The 
hadronic basis of states moving in the "away-side" jet 
will be replaced as before with the sum over all momen- 
tum states of a single antiquark: 



{2tt)32Ei {2tt)32E2 (2n)32k 



E E ■ (59) 



-1S2 



In the last two subsections, the quark attached to 
the electromagnetic vertex is assigned the momentum p, 
while the quark or gluon which materializes from the split 
is assigned the momentum 1. In this subsection, the mo- 
mentum of the fragmenting gluon and quark will be set 
to be 1 and p. The quark attached to the electromag- 
netic vertex will thus have a momentum of p -I- 1. This 



17 



cosmetic reshuffling is solely to ease the extraction of the 
single fragmentation functions. 



matrix elements between the partonic operators ip,Ap 
and the hadronic states will result in the single fragmen- 
tation functions. Absorbing the integrals over undetected 
hadron states and the Fourier integrals, the T matrix el- 
ements may be written as (see Sec. II or Ref. pl|): 



k 




frp iba 



S2-1 



{0\AUy)\S2 - 1,P2)(P2,^2 - l|v4^(0)|0), (61) 

which leads to the definition of the gluon fragmentation 
function at leading twist. The same extraction may also 
be performed for the quark overlap matrix operator: 



m^ix)\s, - i,pi){pi,s2 - i\r,mo), (62) 

which leads to the definition of the quark fragmentation 
function at leading order. Within the coUinear approx- 
imation applied to the above two matrix elements, they 
may be approximated at leading twist as 



9j^p(';p2) 



S'''d^p{l)Tgil;p2), 



FIG. 10: The leading log mixed contribution to the NLO 
modification of the quark dihadron fragmentation function. 

Again, we focus on the leading log portion of the ma- 
trix elements. This essentially restricts our attention to 
the square of the matrix element for the process depicted 
in Fig. ^1 The integrated cross section for this process 
may be expressed as: 



Tq{p;pi) ~ —Tq{p;pi). 



(63) 
(64) 



With the definition of the single fragmentation functions 
and the coUinear approximation, the cross section for in- 
dependent fragmentation may be expressed in a simpli- 
fied form as 



d^Pid^P2 



d'ye''-'^'elg^^2N^C, (60) 



2s 



^^27r<5+((q - 1 - pY)TS:P2)Up-pi)d,.{l) 



y (27r)Hpip2 (27r)8 V T ^^^[(l + P)^^ 2^(l + p)2j^_ 



(65) 



7^(^- W 
(q — k)2 — ie 



x(0|V,(x)bi,5i - l)(5i - l,pi|V;,(0)|0) 



(q - k)2 + ie 



X — {Q\A'',{X)\P2,S2 - 1){S2 - l,P2\Al\{)), 
where two identities of unity 



4; / „ilx 



dH 



This is the leading log contribution to the inclusive cross 
section for the production of two identified hadrons at 
next-to-lcading-order where each hadron emanates from 
the independent fragmentation of a parton. The over- 
lap matrix elements which lead to the definition of the 
fragmentation functions, Tg and Tq have already been 
factored out. This represents a NLO contribution to the 
double fragmentation of the quark emanating from the 
electromagnetic vertex. 

Again, we use a coUinear approximation to isolate the 
leading twist part of the terms inside the curly brackets: 

^I7 —iTT-^'^P'^y^) 



are inserted. As in the preceding subsections the stan- 
dard shift of partonic variables is introduced i.e. k — > 
q — p — I and J d'^k ^ J d^p. The squares of the overlap 



(l + p)2 ' 2 ' (l + p)2 



-Tr 



dpa{l) 



2p+(l-fp)2 ' 2 ' (l + p)2 



(66) 
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After introducing two momentum fractions z'^ and z'2 
through a multiphcative factor of unity, 

1 = dzi d45 (^z; - 8 (^4 - , (67) 

and a rearrangement of the integrals we have the follow- 
ing factorized from of the cross section at leading twist, 



1 ^ e^e\N, 
9 



X Tr 



(4- Y- i^h^^i" 



27r5+((q-l-p)2) 



^2p+(l + p)2' 2 ' (l + p)2 



(27r)4"\-^ /+ 



1=P2/Z2 

P=Pl/2l 
l=P2/22 

P = Pl/^l 



(68) 



It is apparent that the second line of the above equation 
corresponds to the hard cross section of an e"'"e~pair an- 
nihilating via a single virtual photon to a qq pair. The 
third line corresponds to the splitting of the quark into 
a quark and gluon. Note the absence of any 5— function 
maintaining an on-shell condition. This indicates that 
neither the quark nor gluon is being cut. The fourth 
and fifth line indicate the independent fragmentation of 
the quark and gluon into hadrons with the identification 
of a single hadron from each of these sources. The cut- 
vertex structure of this process resembles that of Fig. ^2 
The trace over the Dirac matrix structure of the third 
line may be performed, followed by a contraction of the 
Lorentz indices to obtain: 

z[8{p+/z[+p+/z'2) l + y2 



where the variable y is introduced once again as the in- 
tegral over a (5— function 



dy6\y 



vth- 



(69) 



One notes that the variable y is essentially the ratio of 
the forward light-cone momentum of the offspring quark 
to that of the parent. This leads to the same splitting 
function as that of a quark splitting to a quark and a 
gluon. We may have chosen y to represent the ratio of 
the energy of the gluon to that of the parent quark. This 
would have resulted in a splitting function similar to that 



of the preceding subsection. Changing the order of inte- 
gration between the various ratios y, z'-^^z^-, we define the 
quantities z\ = z'-^jy and Z2 = z^ji^ — y). Note that the 
double differential cross section (Pa/dzidz2 involves the 
ratios zi , Z2 of the hadronic forward light-cone momen- 
tum to that of the parent quark emanating from the elec- 
tromagnetic vertex. Following this we may again switch 
the order of integration. 



[ dz'i I dz2 I dy 
Jo Jo 



1 .1 .1 
dy / dz[ / dz2 

Jo Jo 

\y ^ f'^^'~'^ dZ2 

' Jo y Jo (1 - y) 

dzi / dz2 / (70) 

Jo Jz, y(i - y) 



The integration over the hadronic momenta may now 
be subjected to the same variable transformation as 
demonstrated in Fig. ^ Within the coUincar approxi- 
mation it may be easily demonstrated that 



(1 + P)' 



gjph y(i - y) 

4pfp^ZiZ2 



(71) 



The (5— function introduced in Eq. H69() may also be sim- 
ilarly simplified to obtain. 



Pt/z[ 





"2/(1- 


-y) (vt 






P+ ^ 


-1+\Z2 





y(i -y) \Z2 zi J 

where the quantities p+ and are subjected to the 
condition of constraint introduced in Eq. (|68|l . The 
(5— function is now similar to the second (5— function in 
Eq. (|35| . and may be used to extract the hard cross sec- 
tion erg'' [see Eq. ©]. 

In the last two subsections, there is an integration over 
the transverse momentum of the parent quark emanat- 
ing from the EM vertex. This integration has a coUinear 
divergence that must be absorbed into the renormalized 
fragmentation functions. The light-like null vector n was 
alligned in a direction such that its three-components 
remained opposite to those of ph, the sum of the mo- 
menta of the detected hadrons. By definition n has no 
transverse component. As both hadrons originated from 
the same parton, the preceding condition along with the 
assumption of collinearity of the final hadrons with the 
fragmenting parton, constrained the transverse compo- 
nent of the fragmenting parton to be near vanishing. In 
this case we have again chosen the three components of n 
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the concluding lines of Eq. lO, i.e. 



7 X f , P1+P2 
jr-^O Zl + Z2 

^Ph \ : P 

Pi 






FIG. 11: Next-to-Leading order cut-vertex representation of 
the mixed contribution. 



to be opposite to those of ph- However, the two detected 
hadrons in this subsection emerge from two different par- 
tons which come from a split of the parent quark. As a 
result, whereas both fragmenting partons have a trans- 
verse component {l±,p±) proportional to the transverse 
component of the detected hadrons q± , the parent quark 
is restricted to near vanishing values of a transverse com- 
ponent. In conformity with the previous subsections and 
with the DGLAP evolution of the single fragmentation 
function, the integration over a transverse component 
will remain that of a partonic momentum. The hadronic 
transverse component q± may be trivially related to its 
partonic counterpart thus: 



9-L 



2L 



y{l - y){zi + Z2)' 



(73) 



With these approximations, it may be easily demon- 
strated that the double differential cross section of inde- 
pendent fragmentation assumes the form as sketched in 



dzidz2 



EqqOLg f dl']_ 
'^0 9^ / 72 



Cf 



1 



V \V 



1-Z2 



D„ 



dy 



y(i - y) 

Z2 



1 



y 



(74) 



Where Dg and Dg are the single fragmentation functions 
of a quark and a gluon Q (see also Sect. II). 

The careful reader will have noted that the order of 
integration between y and l± has been switched; whereas 
from Eq. H73|l we see that if q± is independent of y then l± 
is a function of y. There are two ways to resolve this. We 
may ascribe a multiplicative y dependence to the max- 
imum value of q±. Indeed this would be an artificial 
dependence, and would physically mean the inclusion of 
an extra piece of phase space in the l±_ integration. Any 
such multiplicative factor would have no influence on the 
DGLAP evolution equations as these would involve a dif- 
ferentiation with respect to i.e. 



dll 



dl 



-1^- (75) 



Yet another approach may be to view the order of inte- 
grals in Eq. (|74|l as merely symbolic, the differentiation 
with respect to leads to the same evolution equations. 

There remains one last contribution to the inclusive 
cross section for the production of two hadrons from a 
parent quark: this is obtained trivially by switching zi 
and Z2 in Eq. (|74|l . i.e. the hadron with momentum frac- 
tion Zl originates from the fragmentation of the gluon 
rather than the quark. With the addition of the above 
mentioned contribution, we complete the discussion of 
the real radiative corrections to the inclusive cross section 
of same-side dihadron production. The contributions dis- 
cussed in this subsection possess no infrared divergence 
as the y integration is terminated at 1 — Z2 or 1 — zi , con- 
trary to dihadron fragmentation from the single quark 
whose infrared divergence is cancelled by the virtual cor- 
rection. 



V. RENORMALIZED FRAGMENTATION 
FUNCTIONS AND DGLAP EVOLUTION 

So far we have calculated the real radiative corrections 
to the dihadron fragmentation functions. One notices 
that the contribution from quark fragmentation after a 
gluon radiation in Eq. (|52|1 contains an infrared diver- 
gence. Such an infrared divergence will be canceled by 
virtual contributions from interference diagrams such as 
those of Fig. ^] It is well known that in gauge theo- 
ries, in light-cone gauge, the leading log contribution is 
contained solely within the self-energy diagram and not 
shared between the self-energy and the vertex correction. 
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as is the case in the Feynman gauge H- The hadronic 
tensor for such a virtual correction is 



S-2 •' 



{2TTfAEiE2 (27r)4 



27r(5((q-p)2) 



xTr 



Tq{p\Pi,P2){-i^ijp)}^Y{A- i>)i 



, (76) 



where, —i'E^p) respresents the one-loop quark self-energy 
of a quark with four- momentum p. 




cross section is 
d^a 



dzidz2 



^ f"^" dll 



q 

1 1 



y 







(77) 



Note that the transverse integration is over l± and not 
the parent quark momenta. The variable y is defined as 
y = /p^. The overall negative sign should not be a 
cause for alarm, as this is only a part of the total cross 
section. Combining the above equation with Eq. I|52|) 
leads to the cancellation of the infrared singularity as 
?/ — > 1. One can effectively combine the virtual and real 
corrections with a -function, 



q^qg 



Cf 
Cf 



1 



y 



(i-y) 



(78) 



where the 



-function is defined as 



k 



FIG. 12: The leading log self-energy contribution to the NLO 
modification of the quark fragmentation function. 



The hadronic tensor may again be factorized at lead- 
ing twist following much the same procedure as those of 
the last section. There remains the integration over the 
internal gluon momentum 1. The leading behaviour of 
this integral, as in the case of single fragmentation, in 
the part of phase space that includes a pinch singularity 
on the l~ contour, and endpoint singularities on the lj_ 
and Z+ contours. The pinch singlularity may occur only 
in the region where 



0<l+ <p^ 



and 



< ^1 < < 



The derivation of the leading behaviour, which mostly 
mirrors the calculation for single fragmentation func- 
tions, will not be pre sented in full detail here. We refer 
the reader to Ref. |2J|, for details. The final result of the 
self-energy correction to the partial double differential 



(1-2/)+ 



1 - y 



with F(y) being any function that is sufficiently smooth 
a.t y ~ 1. 

In the remaining, we will focus on the non-singlct (NS) 
fragmentation functions for simplicity. In this case, the 
contribution from dihadron gluon fragmentation drops 
out. Summing all three types of contributions from the 
last section, we obtain the NLO contribution to the NS 
dihadron fragmentation function. 



(^1>^2,Q') = D%^l,!'^Z,,Z2) 

rQ" dpi f'dy^Jl + y' 



Pis 

2^ 



2^ 



PI Jz y 

dp, f^-'^ dy 



1 -y y y 



2/(1 - y) 



y 



(80) 



where the leading order fragmentation functions are de- 
fined as matrix elements of field operators in Eq. H37|) 
and (|18() . This has exactly the same structure as we have 
outlined in the parton model in Eq. {T)). Therefore, the 
definition of renormalized fragmentation functions and 
the derivation of the DGLAP evolution can be similarly 
applied here. 

The renormalized dihadron fragmentation function is 
defined as 
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(81) 



We point out, once again, that the scale at which the 
renormaUzed functions are defined is chosen above the 
semihard scale . At this scale, corrections to the renor- 
malized quantities may be evaluated perturbatively. In 
terms of the renormalized dihadron fragmentation func- 
tion, the dihadron fragmentation function at NLO can 
be expressed as 



glog(^)P_*.. 



, (82) 



where we have also used the renormalized form of the 
single fragmentation functions. 

Note that we have introduced a minimum limit to the 
scale fi > in the definition of the renormalized di- 
hadron fragmentation function. If one insisted in choos- 
ing /i < /i^, then the contribution from independent 
quark and giuon fragmentation after a semi-hard split 
would have to be factorized off at the scale fi± . In other 
words the minimum of the second logarithm in Eq. H82() 
would be set at and would only contribute in the 
event that is chosen greater than /i^. In this way, 
we have assumed that only hadron pairs with relative 
transverse momentum q± > fj,± are generated pertur- 
batively from independent fragmentation of two sepa- 
rate partons in the process of a perturbative cascade. 
For two partons whose relative transverse momentum is 
smaller than nonperturbative processes become im- 
portant and their fragmentation cannot be independent 
anymore. We include this part, which also contains a 
coUinear divergence, in the renormalized dihadron frag- 
mentation function. This non-perturbative scale can also 
be considered as the intrinsic relative transverse momen- 
tum of the dihadron fragmentation function and it should 
set the limit of the integration over q± in the matrix el- 
ement definition of the dihadron fragmentation function 
in Eq. (|37|l . If one wants to consider the unintegrated 
(over q±) dihadron fragmentation function, ii± could set 
the initial condition for the qj_ distribution and can be 
used to study the evolution equation of the angular dis- 
tribution inside a jet. For now, this scale will only set a 
limit of the physical scale >> > m1 >> ^qcd 
the DGLAP evolution and will not enter the equation. In 
the ensuing calculation of the evolution of the fragmenta- 
tion functions we will not enter into such subtleties and 
always chose the starting scale /i > 



To include the entire leading log modification, contri- 
butions from all the diagrams outlined in Fig. ^ have 
to be resummed into the scale dependent fragmentation 
functions. These are then differentiated to obtain the 
evolution equation which is given exactly as in Eq. H15(l 
for NS dihadron fragmentation function. 



VI. 



NUMERICAL RESULTS OF NON-SINGLET 
EVOLUTION 



z,=2z. 



log(Q )=0.693 

1.693 
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3.693 
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FIG. 13: Results of the evolution of the non-singlet quark 
dihadron fragmention function Dg^ ''^ (zi , 22 ) , where zi = 2z2, 
from Q2 = 2GeV^ to 109GeV^ See text for details. 

In this section we will study numerically the DGLAP 
evolution of the non-singlet dihadron fragmentation. As 
in many other cases of DGLAP evolution, the solutions 
require an initial condition of the fragmentation functions 
at an initial scale. Such initial conditions, as in the case 
of single fragmentation functions, are non-perturbative 
and are usually constructed from the experimental mea- 
surement of the single inclusive differential cross section 
da/dz, according to Eq. at LO. The evolution of the 
fragmentation functions with the energy scale of the reac- 
tion can then be calculated from the DGLAP equations. 

The absence of any experimental data for two particle 
correlation in e~''e~annihilation forces us to formulate an 
ansatz of the initial condition. We simply use it as a toy 
model to illustrate the DGLAP evolution of the dihadron 
fragmentation functions. We take the LO product of two 
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^2 



FIG. 14: Same as Fig. 1131 except zi — 3z2. FIG. 15: Results of the ratio of the non-singlet quark dihadron 

fragmention function D^^^^ {zi, Z2,Q'^) to the single leading 
fragmentation function Dq^ {zi,Q'^). In this case zi — 2z2 

single fragmentation functions as the initial condition for and = 2GeV^ to lOQGeV^. See text for details. 

the evolution of the fragmentation function, i.e., 



X 0(l-zi-Z2). (83) 

We set the initial condition at = 2GeV^. This cor- 
responds to logQ^ = 0.693. While it may be argued that 
the initial energy is somewhat low for the applicability 
of pQCD in e+e^ annihilation, we consider it as just the 
scale of the momentum transfer while the actual jet en- 
ergy could be sufficiently high. The differential equation 
corresponding to Eq. H15|l is then solved by the simple 
methods of a second order Runge-Kutta numerical esti- 
mation. Results are presented in Figs. ll3l2Ul at intervals 
of AlogQ^ = 1.0 The initial condition is represented by 
the solid black line in all plots. We stop the evolution at 
logQ2 ^ 4 693, which corresponds to ~ lOQGcV^ 

We present results where the leading particle possesses 
a multiple of the momentum fraction of the next-to- 
leading particle. We begin with plots of just the evo- 
lution of the non-singlet dihadron fragmentation func- 
tion at zi = 2z2 in Fig. E| and zi = 3z2 in Fig. ^1 In 
these and all other plots the results are always presented 
as a function of Z2. The results of the evolution are not 
qualitatively different from those of the single fragmenta- 
tion function. Since the sum of the momentum fractions 
are constrained to unity i.e. zi + Z2 < I, we find that 



the fragmentation functions terminate at Z2 ~ 1/(1 -I- r) 
as appropriate. The intial condition, which is merely 
the product of two single fragmentation functions and is 
not subjected to this constraint, does not show this be- 
haviour. It is imposed by hand, in the initial condition, 
that they vanish at and above this value. 

In experiments, one can first identify a hadron as the 
leading hadron inside a jet and use it as a trigger with 
given momentum zi. Then, the associated or the "next- 
leading" hadron distribution inside the same jet corre- 
sponds to the ratio of dihadron and single hadron frag- 
mentation functions, D^^^^'^{zi,Z2)/ D^^{zi). We present 
results for this ratio at zi = 2z2 in Fig. ^1 at zi = 3z2 
in Fig. 1161 and in the extreme cases of zi = Z2 in Fig. 1171 
and zi = 4z2 in Fig. E| It should be pointed out that 
the y-axis in all these plots is linear and not logarithmic. 
Thus, one concludes that the ratio demonstrates little 
qualitative change for a variation of by almost two 
orders of magnitude. In making such an observation, one 
must ignore the difference between the solid line (initial 
condition) and the remaining lines (evolved functions) es- 
pecially at large Z2 . This is due to the fact that the initial 
condition is not subjected to the kinematic constraint as 
Z2 ^ 1/(1 + f)- However, the four plots are visually 
quite different from each other: the ratios of the evolved 
functions display a steady drop as compared to the ini- 
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tial condition as we progress from zi = zi to z\ = 4z2. 
The reader will also note that the maximum of the y-axis 
drops with increasing r. 

As the energy of the e"'"e~annihilation is raised the 
multiplicity must also increase. The cause for this is 
nothing other than the excess energy available for par- 
ticle production. Following the physical picture of frag- 
mentation proposed in Ref . |i22> | , one notes that at higher 
energies it becomes more probable that the two hadrons 
emanate from two causally disconnected sectors of the 
fragmenting jet. If this were the case, then at very 
high energies, the ratio of the double fragmentation func- 
tion to the product of the single fragmentation functions 
{Df'^^ {zi,Z2, g2)/i^Ji (zi, Q2)£,^i (z2, Q2)) should reach 
unity, especially for small values of Zi, zi. This has not 
turned out to be the case as evidenced by the plots of this 
ratio in Fig.^|for z\ = 2z2 and in Fig.[2n|for z\ = 3z2. 
In these plots the ratio deviates from unity at small mo- 
mentum fractions. Whether this is a facet of the choice 
of our initial conditions is as yet unclear. 



VII. DISCUSSIONS AND CONCLUSIONS 

In this paper, we have studied dihadron fragmentation 
functions within the framework of coUinear factorization 
in the high-energy e~ annihilation processes, starting 
with the operator definition. Using the cut-vertex tech- 



nique, we also derived the DGLAP evolution equation for 
the non-singlet dihadron quark fragmentation function. 
We solved the DGLAP evolution equation numerically, 
with a simple ansatz for the initial condition. 

Both the definition in operator formalism and the re- 
sultant DGLAP evolution equations for the dihadron 
fragmentation functions have remarkable similarities to 
the single fragmentation functions. The first type of con- 
tribution to the evolution equations, that from gluon ra- 
diation before the fragmentation of the offspring quark 
(or gluon) into a pair of hadrons, is very similar to the 
corresponding process in single fragmentation functions. 
The second type, unique to the dihadron fragmentation 
functions, comes from independent fragmentation of the 
two offspring partons into two single hadrons of the ob- 
served pair. Since this piece represents the incoherent 
fragmentation of the quark and the gluon, it is well de- 
fined only when the transverse momentum between the 
detected hadrons is large. 

The relative transverse momentum between the two 
hadrons is integrated over in the definition of the di- 
hadron fragmentation function. Hence, we have assumed 
that its non-perturbative contribution, which resides in 
the operator definition, is restricted to an intrinsic trans- 
verse momentum scale, < Hadron pairs with 
Q± > M_L are assumed to be generated only perturba- 
tively. For this assumption to be justified, the semihard 
scale fi± is chosen to be much larger than Xqcd- One has 



24 




' ' ^ ' ^ ' ' 

0.1 0.15 0.2 0.25 

^2 

FIG. 18: Same as Fig. ll5lexceDt zi — 4z2. 



then to assume that the energy scales of the processes in 
question i.e. be much larger than the semihard scale 
/Lt_L, such that a hierarchy of scales, ^qqd m1 
Q^, is satisfied. 

This study is motivated by the observation 0| that the 
same side correlations of two high px hadrons in central 
Au + Au collisions remain approximately unchanged as 
compared with that in p~\-p and d+Au collisions. Specifi- 
cally in this experiment, one measures the distribution, in 
azimuthal angle, of the secondary (or associated) hadron 
T^r~^ with respect to the triggered high px hadron. 
Neglecting the differences in production cross section and 
fragmentation functions for different parton species, the 
integrated yield of the correlation around the peak at 
(j) = should be the ratio of dihadron and single hadron 
fragmentation functions, D'^^''''{zi, Z2,Q^)/D'^^{zi,Q^), 
with zi and Z2 being the momentum fractions of the 
triggered hadron and associated hadrons, respectively. 
To understand the observation in the framework of jet 
quenching, one has to study the medium modification 
to a dihadron fragmentation function due to parton en- 
ergy loss. Since it has been shown in the case of single 
fragmentation functions that medium modification due 
to multiple scattering and induced gluon radiation closely 
resemble that of radiative corrections due to evolution in 
vacuum the DGLAP evolution is expected to yield 
clues regarding the medium modification of the dihadron 
fragmentation functions. 



' ' ^ ' ^ ' ' 

0.1 0.2 0.3 0.4 

^2 

FIG. 19: Results of the ratio of the non-singlet quark di- 
hadron fragmention function D^^'^^{zi, Z2, Q^) to the product 
of the single fragmentation functions Dq^ {zi, Q^)D^^ {z2, Q^). 
In this case zi = 2z2 and = 2GeV^ to 109GeV^. See text 
for details. 



Our numerical results indeed show little change of the 
ratio D'^^'^^{zi,Z2,Q'^)/D'^^{zi,Q'^) as ig varied in 
a wide range of values. The evolution is shown to be 
strongly dependent, however, on the ratio of the momen- 
tum fractions of the two hadrons (r = zxjz^). In the 
results of Ref. Q the ratio r = Z\l z^ is essentially inte- 
grated over all values > 1. In order to relate to the results 
in this paper one essentially must average the effects of 
evolution shown in Figs. (|15I18|) . One will immediately 
note that summing over different values of the ratio z\l z^ 
will lead to the observation of minimal change in the ra- 
tio of the fragmentation functions as a function of the 
of the reaction. 

No doubt, this study is but the first step in this ef- 
fort. In the interest of simplicity, results for the com- 
putationally simpler non-singlet fragmentation functions 
were presented. The results for the evolution of the more 
physically relavant singlet fragmentation functions will 
be presented in a future effort. The DGLAP evolution 
equations for such functions will involve, in addition, the 
splitting of one gluon to two giuons and the coupled dif- 
ferential equations. 

In the above, we have demonstrated the factorization 
of the double differential cross section into a LO hard part 
and a soft piece which encoded the nonperturbativc infor- 
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of this process rests on the assumption that the higher 
order diagrams that have gluon Hnes connecting the out- 
going quark and gluon can be ignored in a leading log 
analysis. The simplest higher order correction of this 
type emanates from the presence of a single gluon line 
connecting the outgoing quark and gluon as illustrated 
in Fig. 

The presence of a Sudakov double logarithm in this di- 
agram signals a leading log contribution to the infrared 
and coUinear divergences of diagrams of the same order. 
Such a contribution will call into question our derivation 
of the second piece of the evolution equation [Eq. (|74|) ]. 
Derivation of this piece of the evolution equations re- 
quired the identification of the leading log contributions 
from infrared and coUinear divergences followed by fac- 
torization and a resummation of these into the single 
fragmention functions. Such a procedure may only be 
carried out if the leading log portions in higher order di- 
agrams are contained solely in the selfenergy corrections 
of the outgoing gluon and quark lines and in real gluon or 
quark emmissions off these lines. A leading log contribu- 
tion from the infrared or coUinear sector of the diagram 
of Fig. 12] will call this procedure into question. In this 
appendix, we demonstrate that in the coUinear sector of 
this diagram, evaluated in light-cone gauge, the Sudakov 
double log is absent. 



FIG. 20: Same as Fig. UHl except zi = 3z2 



mation of converting partons into hadrons. A complete 
proof of factorization requires the extension of the calcu- 
lation to all orders. There also remains the evaluation of 
the medium modifications to the dihadron fragmentation 
bmctions. We will address each of these issues, in turn, 
in future publications. 
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APPENDIX A: SUDAKOV DOUBLE 
LOGARITHMS IN LIGHT-CONE GAUGE 

In Subsect. c of Sect. IV we argued that one of the 
NLO contributions to the dihadron fragmentation func- 
tions would result from the convolution of two single frag- 
mentation functions. The justification for the inclusion 




FIG. 21: A next-to-next-to-leading order correction. 

In this endevour we follow the techniques outlined in 
Ref. for the evaluation of the photon vertex. The 
leading contributions from the diagram of Fig.|^may be 
evaluated as in by the solution of the Landau equa- 
tions. This leads to the presence of two regions of phase 
space where a double logarithm may arise: 
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F - fc'' = aafc^ , k^ = p . k 



(Al) 



In the above, ai,a2 are mere real numbers. These con- 
ditions essentially outline the cases where the gluon A:^ 
achieves collinearity with the outgoing quark with mo- 
mentum or with the outgoing gluon with momentum 
l^^ . Either case produces a double log contribution in the 
Feynman gauge. 

The expression for a quark with momentum p -f 1 split- 
ting into an outgoing onshell quark with momentum p 
and gluon with momentum 1 with the splitting vertex 
corrected by a single gluon with momentum k may be 
expressed as 



u,{p)T;^M,{p + l)eP,{l)^u,{p) 



idaf,{l - k) 



(p + k)2-Kze (l-hk)2+ie 
gP^'ik - 2lf + g°''^{l - 2k)P + g'^P{k + /)" 



k2 



(A2) 



It may be argued (see Ref. Q), that the double log- 
arithm behaviour in vertex diagrams originates on the 
collinear pinch surfaces outlined above in Eq. ljAl|l . The 
behaviour of this above vertex correction in the region 
where the internal gluon momentum k becomes collinear 
to the the outgoing gluon momentum I will now be 
analysed. This corresponds to the second condition in 
Eq. (jAip . The d^fc integration is decomposed into the 
light-cone variables dk^dk~ (Pk±. The focus is on the 
pinch sigularity which results from the two denominators 



1 



1 



k2 



.'(1-k) 



The pinch between the two denominators arises in the 
k^ integration solely in the region < k~^ < l'^. Evalu- 
ating the pole at k^ = k\/2k^ . We obtain a pinch from 
the 

1 



(1 - k)2 + ie 

propagator in the region where k± ■ lj_ = k±lj_, i.e. k 
becomes collinear with 1. To regulate the collinear diver- 
gence we introduce the variables x, Sk±: 



k+ = xl+ 

kj_ = xlj_ + Sk± . 



Evaluating the integrand of Eq. ljA2p at the residue 
of the pole k~ = fc^/2fc+ followed by the substitutions 
outlined in Eq. l|A3p we obtain the split vertex correction 
as 



5kl 



_, , y[,^-vx y)p+i^ „ , 



9a, 



ll{p++l+)^ 

{I - x){lan^+l^no 



(1 -a;)I-n 

""(x - 2)f + g"^{\ - 2x)lP + /"(x + l)r 
dpAl)eT{l). (A4) 



In the above equation all factors of color together with 
multiplicative constants have been absorbed into the fac- 
tor S*". The sum over polarizations of the gluon with 
momentum 1 — k [i.e. the factor dap,{\ — k)] and the 
expression for the glue vertex have been simplified. In 
these and in the rest of the numerator, factors of 5k± 
have been neglected, as the focus is on the region where 
(5fcx — ^ 0. If the numerator, under this approximation, 
turned out to be vanishing, then this would indicate the 
leading contribution to be proportional to 5k±^ and 
as a result no double logarithm and no leading contri- 
bution from this loop would result. Under the relations 
afforded by the polarizations of the light-cone gauge: 



f-d 



F-e^*(0 = 0, 



(A3) 



it may be easily demonstrated that a contraction of the 
Lorcntz indices a, p, (3 lead to the numerator of the r.h.s. 
of Eq. (|A4|) to become vanishing. This result is a property 
of the light-cone gauge. 

It may also be demonstrated, following similar meth- 
ods, that the double logarithm emanating from the re- 
gion of phase space where the gluon k^' becomes collinear 
with the outgoing quark line is also suppressed due the 
vanishing of its coefficient as above. We leave the proof 
of this property to the reader. The above arguments 
demonstrate the vanishing of the leading contributions 
form higher order corrections to the split vertex. This 
property validates our picture of independent fragmen- 
tation. 
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